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We consider the spin-1/2 hamiltonians proposed by Cano and Fendley [J. Cano and P. Fendley, Phys. Rev. 

Lett. 105, 067205 (2010)] which were built to promote the well-known Rokshar-Kivelson (RK) point of quan¬ 
tum dimer models to spin-1/2 wavefunctions. We first show that these models, besides the exact degeneracy of 
RK point, support gapless spinless excitations as well as a spin gap in the thermodynamic limit, signatures of an 
unusual spin liquid. We then extend the original construction to create a continuous family of SU(2) invariant 
spin models that reproduces the phase diagram of the quantum dimer model, and in particular show explicit ev¬ 
idences for existence of columnar and staggered phases. The original models thus appear as multicritical points 
in an extended phase diagram. Our results are based on the use of a combination of numerical exact simulations 
and analytical mapping to effective generalized quantum dimer models. 
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INTRODUCTION 


Construction of low-energy effective models has a long his¬ 
tory in physics, when the original model of interest is too 
complicated to solve or to be treated with perturbation theory. 
This is true for instance in lattice quantum chromodynamics 
or in condensed matter physics. When constructing such low- 
energy effective models, it is important to remember how the 
initial degrees of freedom are transformed, especially if they 
are naturally coupled to experimental probes. 

A canonical example in condensed matter physics is found 
in the context of Mott insulators where the charges are frozen 
and the spins become the relevant degrees of freedom. In that 
case, the electronic Hubbard model at half-filling in the large 
interaction limit can be reduced to the Heisenberg model 


Hj = jJ2 Si.Sj (1) 

(m) 


which only involves spin degrees of freedom. The exchange 
constant J favors antiferromagnetic order in the ground-state 
(GS) when positive and the sum runs over nearest neigh¬ 
bors (NN) spins. Quantum magnets exhibit a large plethora 
of exotic phenomena, especially when the exchange inter¬ 
actions are frustrated. There, magnetic long-range order 
can be prohibited and two NN spins naturally tend to lock 
in a non-magnetic two-sites singlet or valence bond (VB) 

(in) - in»/V2. 

When the physics is local and non-magnetic such as often 
found in frustrated magnets, it has been argued that NNVBs 
can serve as the correct emergent degrees of freedom. It is 
therefore tempting to derive an effective model in terms of 
these short-range objects. However, this is a difficult task 
as VB states are non-orthogonal. Building on this intuition, 
Rokhsar and Kivelson (RK) constructed a model 1 allowing 
quantum dynamics for two-sites short-range objects, which 
can be derived from the first terms of an overlap expan¬ 
sion of the Heisenberg model. The resulting quantum dimer 


model (QDM) 

#QDM = X [-*(!”} (| || + h.C.) + 

«(DCI+II l)(l II) 


( 2 ) 


(here depicted on a plaquette of the square lattice) involves 
(NN) quantum dimers , which are hardcore objects with no in¬ 
ternal structure (in contrast with the antisymmetric VB where 
the spin degrees of freedom can still be probed). Quite im¬ 
portantly, quantum dimer states are taken as orthogonal. An 
important feature of Eq. (2) is the existence of a so-called RK 
point at v — t where the ground-state is exactly known and is 
an equal amplitude state in each ergodicity sector of 77qdm 
(see Appendix A). 

Besides their relative simplicity due to the use of an or¬ 
thogonal basis, there are several practical advantages to use 
QDMs instead of (frustrated) spin models. They are indeed 
more tractable analytically (existence of RK points, 1 some 
QDMs are even exactly solvable 2 ) and numerically (the re¬ 
duced Hilbert space size allows the study of larger systems 
with exact diagonalization and very often no sign problem ap¬ 
pears in quantum Monte Carlo simulations). These techni¬ 
cal aspects only partly explain the popularity of QDMs, the 
other ingredient being the richness of their phase diagrams. 3 
Indeed, several unconventional phase of matters such as Z 2 
or J7(l) liquids or valence bond crystals (VBC) of different 
types (columnar, staggered, plaquette or mixed) have been ex¬ 
hibited in QDMs, 2,4-7 which are hard to find or characterize in 
spin models. 

Quite crucially however, the initial spin degrees of freedom 
are completely forgotten in QDMs as they all are assumed to 
pair in NN singlets. They are yet essential if one wishes to al¬ 
low magnetic order to compete, or just to compare with exper¬ 
imental results which probe two-spin (neutron scattering) or 
four-spin (Raman scattering) correlations. In QDMs, two-spin 
correlations Cij = (Si • Sj) cannot be easily computed while 
four-spin correlators C ijk i = {(Si-Sj)(S k -Si))-Cij.Cki are 
mimicked by (hardcore) dimer-dimer occupancy correlations 
(riijUki) in the sole case where i, j and k, l are NN. 

Clarifying the connection between SU(2) spin models and 
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QDM is a long standing question. One interesting route 
was investigated in Ref. 8 using a decoration scheme or 
more recently using Projected Entangled Pair States (PEPS). 9 
In this spirit, a systematic way of deriving QDMs from 
Heisenberg models was introduced by one of the authors and 
coworkers, 10-12 and the resulting effective models, dubbed 
Generalized Quantum Dimer Models (GQDM), involve terms 
on larger plaquettes. 

In an insightful work, 13 Cano and Fendley (CF) took a 
step in the other direction and introduced two spin-1/2 mod¬ 
els with local spin projectors that display RK-like GS in the 
NNVB basis. In this article, we (i) first investigate the low 
energy properties of the two CF models, and (ii) extend them 
by constructing and studying a continuous family of spin-1/2 
models that closely mimics QDM physics. These spin mod¬ 
els share useful properties with QDMs (exactly known GS, 
good variational approach in a reduced Hilbert space), while 
retaining the true spin degrees of freedom (allowing spin or¬ 
der to set in and computations of spin correlations). Using 
an effective model approach as well as exact diagonalizations 
on both the original spin-1/2 model and its projection on the 
NNVB basis, we show that the obtained phase diagram is in 
correspondence with the one of QDM on the square lattice, 
including explicit signatures of VBC as well of existence of a 
RK-like multicritical point. 

The outline of this paper is as follows. We first briefly recall 
in Sec. I the most important aspects of the CF construction and 
introduce the family of spin models studied along the paper. 
We then consider in Sec. II both CF models, and show using 
exact diagonalizations up to 40 spins and diagonalizations in 
the variational NNVB basis up to 50 spins that they indeed 
possess RK-like ground-states but also have a low-energy ex¬ 
citation structure similar to the RK point of the QDM, in par¬ 
ticular gapless singlet excitations. In the second part of this 
manuscript, we extend in Sec. Ill the CF construction to a con¬ 
tinuous family of spin models and establish a correspondence 
between the SU(2)-invariant spin-1/2 models and the QDM 
phase diagrams for arbitrary v/\t\ values. This is performed 
using numerical simulations as well as analytical arguments 
based on a GQDM mapping of the spin hamiltonian. We fi¬ 
nally close with a discussion in Sec. IIIB 3. Two appendices 
provide details on the exact degeneracy of square-lattice QDM 
models at their RK point (App. A), as well as on the GQDM 
mapping of CF models and their extensions (App. B). 


known as a Klein term, 14 and 


= E^ 3/2 


ryS — /I T)S — 3/2 

4 v 


pS= 3/2 ps=n pS= 3/2 

^ J 7 > ’ 
*4* 


(5) 


where the sums runs over all the square plaquettes p of the lat¬ 
tice and i/f stands for the projection operator on the total spin 
S = fi (with fi = 0 or 1) sector on the cluster c. We consider 
here square lattices with periodic boundary conditions. 

For convenience, we recall that these projectors can be eas¬ 
ily reexpressed in terms of spins operators and thus Eq. (3) 
reinterpreted as multispin interactions : 

£ B..S,)(§ + £ *•■»,) 

0,j)e+ 

^ 1,2 =I(p e m.) 

° = \~ Si-Sj 

1 = 5 + s i .s i . 

The roles of the two terms (4) and (5) are different. The 
Klein term Hi is a projection term designed to annihilate any 
NNVB state. Since this operator is positive, this ensure that all 
NNVB are located in its groundstate manifold. Note that this 
manifold also contains non NNVB states (e.g. VB states with 
aligned VB on the diagonal of plaquettes) - at this stage, the 
ground-state manifold of Hi is therefore largely degenerate. 
This term is gapped to other singlet excitations, and this gap 
is expected to be robust when introducing H^. 

H 2 is designed to introduce dynamics in the GS manifold 
of Hi, and therefore to lift its large degeneracy. The form of 
H 2 will b e such as t0 select the equal-amplitude superposition 
of NNVBs as the ground-state of Hi + Hg. To see this, let us 
detail the dynamics introduced by H 2 on a gi yen NNVB state 

W : 



I. CANO-FENDLEY MODELS AND THEIR 
GENERALIZATION 

We recall here the main ingredients of the Cano-Fendley 13 
construction to obtain spin-1/2 models with ground-state(s) 
that are equal-amplitude superpositions of NNVB coverings, 
here on the square lattice. The Cano-Fendley models, 

H M = Hi + H 2 M = 0,1 (3) 

contain two terms: 

Hi=J2 P + =5/2 ’ ( 4 ) 

+ 


1. Case 1: p contains 0 or 1 dimer in \^) (non-flippable 
plaquette). This implies that at least one dimer is hosted 
by a corner sharing neighboring plaquette. In that case 


either p 5 ,- 3 / 2 , pS- 3 / 2 , pf- 3 / 2 or ps-3/2 


'4- m - W 




annihilates 


the state. 


2. Case 2: p contains 2 dimers (flippable plaquette). We 
note that 

<6) 

where Ii j 3 is the permutation operator on sites i and j. 
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It is easy to show that: 



<VOHG> 

(7) 


n^OHQ} 

(8) 


SO) -ID) 

(9) 


SJOHO)- 

(10) 

Hence, 



pS=fi _|_ pS=fi 
. ^ , 

)lO)=< 2 »*- 1 >lu]) + O) 

(ID 


)luH 2 '*- 1 >IO)+ |p) 

• (12) 


At this stage, it seems from Eq. (11)-(12) that 7^ perfectly 
mimicks the QDM Eq. (2) in the variational NNVB basis with 
v = 1 and t = 1 (respectively t = — 1 ) for /i = 0 (respec¬ 
tively /I = 1). However, it is important to remark that the 
states involved in the r.h.s. of Eq. (11)-(12) are not invariant 
under the application of P , or P , ' and these oper- 

ators induce further reconfigurations of NNVB states outside 
the flippable plaquette. This will lead to subtle differences be¬ 
tween the two models that will be discussed in the following, 
together with the precise correspondence between CF models 
and their QDM counterparts. 

Note as well that on the square lattice the sign of t is irrel¬ 
evant for the QDM as it can be gauged away. Therefore, both 
CF models are designed to have the equal-amplitude states 
in the NNVB non-orthogonal basis 15-17 as zero-energy GS. 
Similarly to the situation of the RK-point of the QDM, all 
equal-amplitude states in each ergodicity sector are valid GS, 
therefore the minimal degeneracy of the CF models is equal 
to the number of ergodicity / topological sectors for NNVB / 
dimer coverings of the square lattice (see Appendix A). Note 
that we used the term ’minimal’ as other GS may arise due to 
the Klein term which also annihilates non-NNVBs states on 
a square lattice with periodic boundary conditions. In Sec. II 
we present numerical evidence, using exact diagonalizations 
both in the full S z basis (allowing to study magnetic and non¬ 
magnetic states) and in the variational NNVB basis (expected 
to be relevant there), that the models proposed by Cano and 
Fendley are indeed realizations of SU(2)-invariant spin-1/2 
models that not only share the GS properties but also the low- 
energy structure with the QDM at the RK point. 

In Sec. Ill of this article, we investigate, using various nu¬ 
merical and analytical methods, a generalization of Cano and 
Fendley proposal allowing to mix both CF models H° and 

h 1 : 


n K ’ e = KWi + -F (cos e U° 2 + sin 6 Hi), (13) 
n[6) V / 

where the normalization factor n(6) = |cos#| + |sin0| is 
introduced for convenience. The two CF models are recovered 
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FIG. 1. (color online) Non-NNVB zero energy singlet (resp. zero 
energy triplet) eigenstate of the Klein term: the central red bond is to 
be interpreted as a singlet (resp. a triplet). 


as n° = 7^=1,0=° and i0 = t^=M=P 2 . Note that ^K.e 
also conserves the SU(2) spin symmetry of the CF model, as 
well as square-lattice spatial symmetries. 

The motivation to introduce H K}6 is to see whether this 
family of spin-1/2 models can host phases similar to those ob¬ 
served in the phase diagram of the conceptually simpler QDM 
on the square lattice (Eq. 2), as the mixing angle 6 is tuned. 


II. CANO-FENDLEY MODELS: 0 = 0 AND 0 = tt/2 
A. Prefactor of the Klein term 

Before going into the detailed numerical study of 'H K ’° and 
j s important to investigate the influence of the pref¬ 
actor K of the Klein term Hi in Eq. (13) as it plays an im¬ 
portant role to isolate the NNVB states from other singlet or 
triplet states in the low-energy manifold. 

First let us consider the Klein term alone. On a square 
lattice with periodic boundary conditions, the Klein term is 
known to be not perfect 18 in isolating the NNVB manifold, as 
other singlet states and even triplets also have zero energy 7 . 
Indeed, it is quite easy to verify that the states depicted in 
Fig. 1 are zero energy non-NNVB singlet or triplet eigenstates 
of the Klein term, as the total spin of any group of five spins 
formed by a site and its four nearest neighbors is lower than 
3/2. 

Since the Klein term has no spin gap, it is important to show 
that there exists a regime of value for K in the presence of 
H% for which (i) a sizable spin gap opens in the thermody¬ 
namic limit and (ii) the other interactions included in H% al¬ 
low NNVB configurations to gain further energy by resonat¬ 
ing, and therefore to separate from other zero-energy states of 
the Klein term. 

In the original CF proposal, its amplitude was set to 1, but 
we find that it is necessary to consider a larger value. We show 
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FIG. 2. (color online) Spin gap versus amplitude K of the Klein term 
for the 'H k, ° model on different square clusters with N sites. Inset: 
scaling of the spin gap vs 1/N for two values of K. 


in Fig. 2 the evolution of the spin gap vs the amplitude K 
for the model i-L K ’ e=0 . The finite-size behavior is not always 
monotonous, but a finite-size scaling analysis (see inset) indi¬ 
cates a finite spin gap in the thermodynamic limit quite con¬ 
vincingly when K > 10, while data with K = 1 are more am¬ 
biguous. Since the iterative exact diagonalization techniques 
that we use for the largest samples become slower when K is 
too large, we find it practical to fix its value to K = 10, which 
will be enough to ensure that all models studied further have 
a finite spin gap. This value K — 10 will be implicitly used in 
the rest of this work. We will show in the next section that it 
is also large enough to eliminate non-NNVB intruder singlet 
states from the low energy manifold. 

As a final remark, let mention that in the NNVB approach 
(consisting in a ED of the spin hamiltonian in the subspace 
spanned by NNVB states) that we use in conjunction with ED 
in the S z basis, the two points mentioned above are naturally 
fulfilled since, by construction, neither triplets nor non-NNVB 
states enter the computation. In a sense using this approach is 
equivalent to forcing the Klein term to be perfect. 


B. Ground-state degeneracy and properties of CF models 

Here we discuss the degeneracy of the ground-states for the 
U K $ and models. All ground-states at these points 

have an energy which is exactly 0. Since one can build an 
equal-amplitude ground-state in each ergodicity sector, we 
expect at least the same number of ground-states as at the 
RK point of the QDM. In our exact diagonalization compu¬ 
tations, two cases need to be distinguished for the computa¬ 
tion of the degeneracy of these zero-energy states: (i) we can 
compute exactly the full spectrum for diagonalizations in the 
variational NNVB basis as well as in the S z basis for small 
enough samples (N < 20), (ii) when we cannot compute the 
full spectrum in the S z basis even when decomposing it into 
symmetry sectors, we use Davidson algorithm to compute ex¬ 
actly the (possibly degenerate) lowest-lying states. 


N 16 20 26 32 36 40 50 
GS deg. 17 13 16 69 41 29 47 


TABLE I. GS degeneracies of the QDM at the RK point. Except for 
N = 16 where spurious states enter the GS manifold, these degen¬ 
eracies are identical for the CF models 'H K, ° and 'H K,7V ^ 2 . 


We find that the number of zero-energy GS of the CF and 
QDM models (at the RK point) are exactly the same (see table 
I), except for the sample N = 16 which exhibits 4 spurious 
ground-states (3 singlets and even a triplet) leading to a de¬ 
generacy of 23 instead of 17. 

The NNVB calculation does not take into account (by con¬ 
struction) triplets or long range singlets and does not display 
the TV = 16 anomaly in the GS degeneracy. This allows us 
to interpret spurious singlets as a small size artifact: on such a 
small cluster, long-range singlets cannot exceed 4 lattice spac- 
ings and are still sensitive to the 8-spin CF terms, but this 
anomaly does not persist for larger clusters. The degenera¬ 
cies of both 'H k ’° and fi K ^/ 2 in NNVB calculations are thus 
matching the ones of the QDM at RK point for all cluster sizes 
from N = 16 to N = 50. We also explicitly checked that each 
GS is indeed of equal-amplitude type in the NNVB basis, ex¬ 
actly as the QDM GS in the orthogonal quantum dimer basis. 
The spatial quantum numbers are also identical. 

The physical properties of the GS of the CF model and the 
QDM are slightly different than the ones of the QDM. First of 
all, spin correlations can be computed (opposite to the QDM 
as discussed in the introduction) and are known since a long 
time 19 to be exponentially decaying Cij oc exp(— \i — j |/£) 
with £ = 1.35. 16,17 More recently, it was discovered that four- 
spin correlators Cijki (with z, j and l NN) rather decay al¬ 
gebraically with an exponent unnvb — 1.16, 16,17 i.e. more 
slowly than in the QDM case where «qdm = 2 for the decay 
of dimer-dimer correlators. 20,21 These critical four-spin cor¬ 
relations in the ground-state strongly suggest 13,16 that the CF 
models are gapless with respect to singlet excitations, while 
remaining gapfull for the spin excitations (at least for K > 10, 
see Fig. 2), forming an unusual spin liquid. This is what we 
investigate in the next section by looking at the low-energy 
singlet spectra. 


C. Spectrum of CF models 

Spin gap. A shown in Fig. 3 (left and middle panels), for 
both 'H k, ° and 'H K ’ 7r / 2 (with from now on, K = 10), the 
triplet energy scale is well decoupled from low energy singlet 
excitations. For 7i K, °, our data are compatible with the ex¬ 
istence of a sizable spin gap in the thermodynamic limit. Its 
precise value is not an important property of the construction 
since, as shown in the previous subsection, it can be easily 
modified by an appropriate tuning of AT. A similar behavior is 
observed for i~i K,7r / 2 despite a less clear finite size scaling. 

Singlet excitations. In order to characterize singlets exci¬ 
tations of i~L K ^ and 7i K,7r / 2 , we performed momentum re¬ 
solved exact diagonalizations up to N = 36 sites and NNVB 
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FIG. 3. (color online) Triplet and singlet gaps on various clusters for 
ft K, ° and 7i K,7r / 2 , using K — 10. For comparison, we also plot the 
low-energy excitations of the QDM on the square lattice. 


diagonalizations up to N = 50 sites. Before entering into the 
discussion of the results, let us make an important comment 
on the NNVB approach in the context of the two models stud¬ 
ied here. The NNVB calculation, as it relies on a truncation 
of the singlet subspace, leads to a variational approximation 
of the GS energy in each symmetry sector. The singlet gap 
being a difference of two variational energies, its value cannot 
be strictly interpreted neither as an upper nor as a lower bound 
in general. Here, the situation is different as the zero-energy 
ground states are exactly captured by the NNVB approach. 
The singlet gap value obtained from this approach is there¬ 
fore also variational, providing an upper-bound for the exact 
singlet gap. 

The lowest singlet excitation, represented on Fig. 3 (left 
and middle panels), shows for both models a non-monotonous 
scaling as the system size increases. This can be ascribed to 
the fact that the first singlet excitations is not located at the 
same k-point from one cluster size to another. For each size 
where ED and NNVB can be directly compared, both calcula¬ 
tions leads to the same k-point location of the gap and the data 
show that the NNVB approximation provides a very good up¬ 
per bound of the singlet gap. Nevertheless, even by perform¬ 
ing NNVB diagonalizations up to N = 50 sites, the finite size 
scaling is still rather inconclusive although compatible with 0 
in the thermodynamic limit. 

In contrast, the situation for the (7r, 0) singlet gap is much 
clearer : both ED and NNVB data exhibit a linear scaling in 
1/N which convincingly extrapolates to 0 for infinite systems. 
H K, ° and are thus gapless in the singlet sector. 

QDM comparison. The low energy singlet spectrum struc¬ 
ture discussed above, is very reminiscent of the QDM spec¬ 
trum at the RK point 5 and in particular the presence of the 
piOns (i.e. ( 7 r, 0)) low-lying branch. The right panel of Fig. 3 
displays the finite-size scaling of the lowest gap and the (7r, 0) 
gap for the QDM up to N = 80 sites. This indeed confirms 
that, at a qualitative level, f{ K ^/ 2 and the QDM at the 

RK point share similar low-energy excitations, where even the 
finite size scaling non-monotonicity is reproduced. In fact, be¬ 


cause of the quadratic dispersion of the resonons at ( 7 r, 7 r), one 
expects a gapless spectrum at any momentum k. 22 

Comparison of energy scales in and H K,7r / 2 . Even if 
the two models mimic the same dynamics for nearest neighbor 
SU(2) dimers (see section I and in particular Eq. (11)-(12)), 
the two models remain different due to the presence of spin- 
3/2 projectors on the opposite corners of square plaquettes 
that involve non-trivial NNVB reconfigurations outside flip- 
pable plaquettes. Comparing H K,0 and fi K ^/ 2 data reveals 
(see Fig. 3 left and middle panels) that the main difference is 
the energy scale at which singlet excitations occur, the gaps 
of 'H k, ° being an order of magnitude smaller than those of 
f-^K,ir/ 2 ' intriguing difference, that seems difficult to ex¬ 
plain with simple arguments, will be addressed in the next 
part using a GQDM mapping. We will show that it can be as¬ 
cribed to the presence of effective larger loops processes cor¬ 
rections to the QDM dynamics which are implicitly encoded 
in the spin projectors and are different in magnitude for the 
two models. 

Conclusion. Our results indicate that 7i K, °, i~i K,7r / 2 and the 
QDM at the RK point have a rather similar low-energy struc¬ 
ture. It is therefore tempting to extend the CF construction to 
build a family of SU(2) invariant spin-1/2 models which low- 
energy properties will mimic those of the QDM at an arbitrary 
value of v/\t\. 


III. ARBITRARY 0 MODELS 

Having established the physical properties of the two mod¬ 
els f[ K ^/ 2 that were introduced by Cano and Fend- 

ley, 13 we now move to arbitrary mixing angle 0 realizations 
and show that they can host any phase realized in their QDM 
counterparts. To motivate this relationship, we first derive an 
effective dimer model using the GQDM formalism. We then 
perform direct numerical simulations of 'H K,e for three differ¬ 
ent values of 6. As discussed in Sec. II A, we need for these 
simulations to consider a proper amplitude K for the Klein 
term, and we again found that setting K = 10 is sufficient (i) 
to ensure a finite spin gap for all 0 that we will consider in this 
section (data not shown), and (ii) to be handled numerically 
without introducing a too large cost. 


A. GQDM approach 

GQDM mapping. In some recent works, 10-12,23 a general 
scheme allowing to perform analytically the projection of any 
SU(2) spin model on the non-orthogonal NNVB manifold was 
developped and applied to various cases. The effective models 
obtained are defined in the orthogonal basis of nearest neigh¬ 
bors hardcore dimers and only involve (i) potential (diagonal) 
terms sensitive to the possibility of flipping dimers along a 
plaquette of a given shape (ii) kinetic (off-diagonal) terms that 
actually perform the flip along a plaquette. This type of gen¬ 
eralized QDM hamiltonian differs from the QDM Eq. (2) by 
the fact that flips are no longer limited to square plaquettes. 
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The shape of the considered plaquettes depends on the lat¬ 
tice and on the precision required for the projection. Since 
the amplitude of plaquette terms in the effective hamiltonian 
decreases exponentially with the number of dimers involved 
in the flipping process, it is generally sufficient to truncate the 
hamiltonian retaining only the smallest loops contributions. 
Recently, other schemes to derive GQDM from Heisenberg 
model have been proposed claiming to refine the procedure 
in the context of the kagome and square-kagome antiferro- 
magnets . 24,25 One has still however to perform a direct nu¬ 
merical comparison with the parent spin hamiltonian to un¬ 
derstand whether its physical properties are well reproduced 
by the effective models. Here we will use the original GQDM 
scheme 12 and make a precise connection between the spin 
hamiltonian 'H K ' 6 Eq. (13) and the QDM Eq. (2) that is fully 
confirmed by exact diagonalizations. 

On the square lattice, the smallest plaquette terms are : 

I IXZI + IZXI II- (14) 

I II 11+ IDO (15) 

Here, we will also consider next-leading plaquette terms : 


□ 


□ 


□ 


I 




ZXI I 11+II I 
ZXI Zl+ll I 


IXZ I 
IXI z 
IXI I I 


(16) 

(17) 

(18) 
(19) 


The effective GQDM hamiltonian is defined as 
q-i/ 2 'h k , 0 q-i/ 2 ^ w j iere q an( j are the overlap and 
hamiltonian operators : 


0 = a 2 

□ + a 4 ' |~M 

+ larger loops , ( 20 ) 

n K ’ e = h n a 2 p 

+ 

□ 

□ 

+ hcja 4 □ ( 21 ) 

+ h 1 — 

] □ 

~h hi.jjQ: 2 


[-□ + h\ m j—□ 


+ larger loops , 


with a 2 = 1/2 (this value originating from the non orthogo¬ 
nality of the NNVB basis and the overlap rule), and the values 
of , h \... n are determined by 'H K,e . It is important to note 

that in the GQDM scheme , 12 the above expressions are not se¬ 
ries in a 2 but expansions on GQDM operator basis. In other 
words, no assumption is made on a. 

The final effective Hamiltonian for EL K,e takes the form 


where the analytical expressions for parameters (which all 
depend on the mixing angle 0) have been reported to the Ap¬ 
pendix B (see Eqs. B10-B15) for conciseness. A few specific 
values are discussed in the following (see Table II). 

Some remarks on i~i K,0 and i~i K,7r / 2 . Interestingly enough, 
the two models introduced by Cano and Fendley , 13 which cor¬ 
respond to 0 = 0 and 7 r/ 2 , exactly map to the same RK hamil¬ 
tonian 1 1 4 1 = U 4 = 1/4 if truncated to the first relevant terms 
£4 and t> 4 . Pushing the mapping further reveals differences 
between H K,0 and H K,7T ^ 2 , as they for instance introduce dif¬ 
ferent dynamics (such as the t$ term) in the NNVB manifold. 

As to prove the ability of the GQDM mapping to capture 
the properties of EL K,e , we first investigate whether the cor¬ 
rections £6($), ^6(0), £ 4 ( 0 ), ^ 4 ( 0 ) to the QDM can explain the 
difference in the low-lying singlet excitations energy scale be¬ 
tween 'H k, ° and 'H k ’ 7T / 2 observed in section IIC. On Fig. 4 
the ED spectra of GQDM models corresponding to and 
7 i K ^/2 are represented for a 36-site cluster in each symme¬ 
try sector. The data show indeed a clear difference in energy 
scales of the lowest excitations, very similar to the results ob¬ 
tained with ED and NNVB calculations on the original spin 
models (see Fig. 3). 

Exploring v/\t\ by varying 9. Using the expressions of 
v±(6) and £ 4 ( 0 ) reported in Eqs. (BIO) and (Bll) of Ap¬ 
pendix B, we represent on Fig. 5 the variation of the ratio 
7 ^ 4 /1£ 4 1 as a function of 0. As can be seen from this figure, 
any value of the ratio ^4 /1 £ 41 can be reached by an appro¬ 
priate tuning of 6. The value of ^4 /1 £ 4 1, together with the 
knowledge of the QDM phase diagram, serves as a guideline 
to understand which phases are accessible for the ground-state 
of 'H Kjd for a given 0. The phase diagram of the square lattice 
QDM exhibits VBC states 26-28 such as columnar for v/\t\ < 0 
and staggered for v/\t\ > 1. In the intermediate region, a 
plaquette phase has been suggested to exist , 26-28 or a mixed 
columnar-plaquette , 29 but its presence has recently been ques¬ 
tioned and the columnar phase seems to extend up to the RK 
point . 30,31 

Of course larger loops processes, in particular £g(0), must 
also be taken into account and can somehow modify this sim¬ 
ple picture. A careful numerical study of H K,e should there¬ 
fore be pursued, which will be done using both exact and 
NNVB diagonalization, together with an exact diagonaliza- 
tion study of the GQDM counterpart. In the following we will 
concentrate on three particular values of 0 denoted 6a, 9b 
and 9c (see Table II and Fig. 5), for which we can expect re¬ 
spectively columnar VBC, putative mixed plaquette-columnar 
VBC and staggered VBC phases from the ^ 4 /1£ 4 1 QDM proxy. 


B. Low-energy excitations and dimer correlations 


U eG ’ 9 = -U{9) 
+ v 6 (0) 


+ V4,(0) 

-m 


■ * 6 ( 0 ) 


FT 


( 22 ) 


!-□ +v' 4 (0) 


+ larger loops 


In order to characterize possible symmetry breaking, it is 
useful to investigate the low-energy levels quantum numbers 
on a finite cluster. The VBC phases we expect for models A, 
B and C will have clear signature in the low-energy spectra as 
singlet excitations with definite quantum numbers that should 
collapse onto the ground-state in the thermodynamic limit . 29 











































7 



11 * 

0 

h 1 

6 — k /2 

A 

= — arccos (-^) 

B 

9b = - arccos (^) 

C 

9c = arccos (^) 

u 

0.25 

-0.25 

-0.125 

0.25 

-0.125 

V4 

0.25 

0.25 

-0.25 

0.125 

0.25 


-0.0885651 

0.202914 

0.015695 

-0.117152 

0.130044 

Vq 

-0.031203 

-0.0308274 

0.0311091 

-0.0156954 

-0.0309213 

*4 

0.063808 

0.0742716 

-0.0664239 

0.0292881 

0.0716557 

U4 

0.0364349 

0.0779137 

-0.0468046 

0.00784772 

0.067544 

V4/\u\ 

1 

1 

-2. 

0.5 

2. 


TABLE II. GQDM amplitudes (see text) for several specific values of the mixing angle 6. 




FIG. 4. (color online) Exact spectra of the GQDM models at 6 = 0 
(top panel) and 6 — 7 t/2 (bottom panel) for a 36-site cluster. States 
are labelled according to their momentum (one data set per nonequiv¬ 
alent k point). For k = (0,0) (respectively (7r, 7 r)), we also add a 
label corresponding to the irreducible representation of (respec¬ 
tively C 2 V ) point group symmetry. 


Moreover, dimer ordering can also be checked directly by 
computing (connected) dimer-dimer correlations on the finite- 
size ground-state. We describe in detail below the low-energy 
physics of model A, B and C, based on the collection of our 
numerical data on a N = 36 square cluster presented in Fig. 6 . 



FIG. 5. (color online) Variation of U 4 /IGI as a function of 6. Five 
specific values of 6 are displayed as red dots corresponding to 'H K, °, 
'H k,7r / 2 and 3 other values of ^ 4 /1^41 studied in the main text. 


1 . Model A using 6 a 

According to table II, such a combination would map onto 
a related QDM model with ^4 /1 ^ 4 1 = —2 which is known 
to be deep in the columnar phase . 1 Indeed, Fig. 6 (a) shows 
the low-energy spectra obtained on a N = 36 cluster by ED 
and NNVB simulations for the spin model as well as for the 
GQDM for the point 6a • In all spectra, one can easily identify 
four almost-degenerate ground-states with quantum numbers 
corresponding to columnar order (i.e. (0, 0) A 1 , (0, 0) B\ and 
two-fold degenerate ( 7 r, 0 ) Bi) and a sizable gap above them. 
It is also remarkable that a large part of the low-lying singlet 
excitations structure is common between the ED and NNVB 
computations and is also reproduced by the GQDM. 

The top panel of Fig. 6 (a) represents graphically the dimer 
correlations C ijk i = ((S* • Sj)(S k • S;)) - (Sj • Sj)(S k ■ S t ), 
taken with respect to a reference bond with NN sites (i. j) 
denoted with black dots, and where the position of the NN 
sites (&, /) is varied. This calculation is done using the ground- 
state in the NNVB basis. One clearly observes a columnar 
pattern, in agreement with the analysis above. 
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FIG. 6. (color online) Phase diagram of the QDM (second top panel), dimer correlations (first top panels, for clarity only NNYB data on 
N = 36 cluster) and low-energy spectra of 7~L K,d from exact (fourth panels from top) and NNVB diagonalization (third panels from top), and 
of the related GQDM (fifth panels from top). For dimer correlations, blue (red) bonds denote positive (negative) correlations, while the width 
of each bond is proportional to Cijki (see main text for definition). The energy spectra are represented as a function of the wave-vector and 
point-group symmetry groups of the N = 36 square lattice cluster (see Fig. 4 for labels). The squares surrounding the wave-vector and point 
group labels highlight the sectors for which low-lying states are expected in order to form a columnar state (blue squares of the left panel), 
staggered state (blue squares on the right panels) and a columnar or plaquette state (respectively blue and red squares on the middle panel). 
The corresponding energy levels are pointed out with a circle. ED data are also labeled differently for S z — 0 and S z — 1 (black circle and 
red square respectively). 
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2. Model C using Oc 

One expects to engineer this way a spin model mimicking 
the QDM with ^ 4 /1 ^ 4 1 = 2 which is known to possess an 
exact staggered ground-state. Let us mention that other non 
flippable zero-energy states can be constructed so that the ac¬ 
tual degeneracy of the QDM is extensive: for instance, we 
find 28 zero-energy states on N = 36 cluster for the QDM 
model with these couplings. Those include not only the 4 
states (0, 0 )Al, (0, 0)L>i, ( 7 r, tt)Ei and ( 7 r, tt)E 2 forming the 
staggered ’’pure” states but a total of 4(2 L / 2 — 1) states (with 
N =s L x L) all belonging to the the diagonals of the Bril- 
louin zone (see Appendix A). The overall degeneracy of 28 is 
decomposed in the following way: 6 states with momentum 
(0,0), 8 in (±7 t/3, ±tt/3), 8 in (±2 tt/3, 2 ± tt/ 3) and 6 in 

(7T,7r). 

In Fig. 6(c), one can readily observe that the low-energy 
spectra of the spin-1 / 2 model (both with ED and in the NNVB 
basis) as well as of the GQDM perfectly display this exact 
same degeneracy. Again, the first non-zero singlet excitations 
are also identical between ED and NNVB, and well repro¬ 
duced by the GQDM. 

With such a ground-state degeneracy, the correlations de¬ 
pend on the chosen ground-state. Taking the fully symmetric 
ground-state (i.e. with momentum (0, 0) and fully symmetric 
with respect to the C± v point group), the existence of a stag¬ 
gered phase is confirmed by a direct computation of the dimer 
correlations in Fig. 6(c). 


3. Model B using 0 b 

As a last example, we studied H K ^ B (see table II) which 
is expected to mimic QDM at v±/\t ^1 = 1/2, i.e. in a region 
where the columnar order is not fully established yet even us¬ 
ing large-scale quantum Monte-Carlo simulations. 30,31 

Fig. 6(b) presents the low-energy spectra where we observe 
in ED and NNVB computations and GQDM spectra, that both 
set of states expected to construct columnar or plaquette states 
in the thermodynamic limit are present at low-energy. Note 
that other states in different sectors not expected to contribute 
to columnar or plaquette ordering (such as e.g. k = ( 77 /3,0) 
or k = (27 t/ 3,0)) also have low energy. The competition 
on finite clusters between these different states make for a 
slightly not as good agreement between the three types of 
spectra for this model with 0b, as compared to 0a and 0c 
where the low-energy picture is clearer. 

A similar observation of low-energy states in competition 
(between the columnar and plaquette states) can be made on 
the related QDM model and had led to the conclusion that it 
could be a signature of a mixed phase. 29 However, in the case 
of the QDM, there is no clear gap above these eight states 
and refined ED and QMC studies have found no evidence for 
such mixed phase: 30,31 it is instead indicative of a U(l) regime 
where no crystalline order is frozen yet (on the available sys¬ 
tem sizes). Presumably, the system sizes that can be reached 
are too small to give any reliable information on the thermo¬ 
dynamic limit because of the existing diverging length close 


to RK point. 30,31 

While in principle, the supplementary terms ( t §, vq etc) in¬ 
duced by the non-trivial singlet dynamics in H K,e could fa¬ 
vor one particular type of ordering, we cannot make either 
any definitive statement on the behavior in the thermodynamic 
limit for the point 0b, based on our finite cluster simulations. 
This is confirmed by computations of correlations in Fig. 6(b), 
which indicate also a kind of mixed phase with both columnar 
and plaquette patterns (see Ref. 32 for a related discussion). 


CONCLUSION 

Using numerically exact simulations as well as an ana¬ 
lytical mapping to a generalized quantum dimer model, we 
have studied in detail the low-energy properties of fully-SU(2) 
symmetric spin-1/2 models introduced by Cano and Fend- 
ley on the square lattice. 13 We confirm that the original mod¬ 
els in Ref. 13 indeed admit equal-amplitude superpositions of 
nearest-neighbor valence bond states, which is what they were 
designed for. Besides the corresponding ground-state degen¬ 
eracy, our calculations show that the low-energy physics is 
also in correspondence with the one encountered in the quan¬ 
tum dimer model. More precisely, the CF models admit gap¬ 
less spinless excitations while having a spin gap in the ther¬ 
modynamic limit (for sufficiently large Klein term), thereby 
forming an unusual type of spin liquid (the ground-state dis¬ 
plays power-law dimer correlations). 

Mixing the two models, provides a SU(2) spin-1/2 Hamil¬ 
tonian analog to QDM : we find that we can reproduce 
quite trustfully the phases encountered in the quantum dimer 
model, namely columnar, staggered (and apparent plaquette or 
mixed) valence bond crystals in these spin-1/2 models. This 
strongly suggests that the CF models are multicritical points 
similar to the Rokhsar-Kivelson point of the quantum dimer 
model, and that the unusual spin liquid encountered there is 
not generic. 

Our work can possibly be extended in several different di¬ 
rections. First, it would be interesting to see whether tech¬ 
niques specific to the RK point of the QDM, could also be 
used for CF models: for instance, we could try to construct 
variational states using a single mode approximation to un¬ 
derstand more precisely the spinless excitations that become 
gapless in the thermodynamic limit. 5 Also, it would be inter¬ 
esting to see whether exact dynamics could be performed, at 
least in the NNVB sector, using the techniques introduced in 
Ref. 33. 

Secondly and away from the RK point, it would be inter¬ 
esting to perturb the columnar crystal (as e.g. found for 0a) 
with an antiferromagnetic Heisenberg exchange between NN 
spins, which is now possible since we directly have access to 
spins. The possibility of a deconfined quantum critical point 34 
between the Neel and VBC phases can be explicitly studied, 
albeit on the relatively small sizes available to exact diagonal- 
ization. 

Finally, a construction similar to the one introduced by CF 
for the square lattice, can be performed on non-bipartite lat¬ 
tices, where the QDM exhibits Z 2 liquid phases close and 
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W x 


FIG. 8. (color online) Example of the sign convention used to deter¬ 
mine the winding numbers (W x , W y ) of a dimer configuration (see 
text). 


FIG. 7. (color online) In blue, one staggered dimer configuration 
with winding (0, L/2). In red, by flipping dimers along two diag¬ 
onal lines, one can generate another non-flippable configuration in 
the sector (2, L/2 — 2). Inset shows the allowed winding sectors 
inside the 0 < \W X \ + \W y \ < L/2 square, as well as the sectors 
corresponding to these two configurations. 

at the RK point . 2,4 For instance, we are currently checking 
whether the ground-state of a CF-like model on the kagome 
lattice is adiabatically connected to the one of the Heisenberg 
model, which has been recently argued to host a Z 2 spin liq¬ 
uid . 35-37 
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Appendix A: Ground-state degeneracies of the RK point on the 
square lattice 

Let us consider the QDM at the RK point on the square lat¬ 
tice. It is well known that the equal superposition of all dimer 
configurations is a ground-state of the hamiltonian. However, 
due to the existence of several ergodic sectors, there are in 
fact many ground-states and the purpose of this appendix is 
to compute this degeneracy explicitly, something was never 


presented to the best of our knowledge. 

For a given LxL square lattice with periodic boundary con¬ 
ditions, dimer configurations can be labelled using their wind¬ 
ing numbers (W x , W y ), 0 < \W X \ + \W y \ < L/ 2 , see Fig. 7. 
W x (respectively W y ) can be easily computed by summing, 
on each link (with an alternating sign) cutting a line along x 
(respectively y) direction, the contribution ±1 if the bond is 
occupied or 0 if it is empty (see Fig. 8 for the sign conven¬ 
tion). It is straightforward to see that these labels are invariant 
with respect to any local rearrangement of the dimer, i.e. they 
are conserved quantities for the Hamiltonian. So, there are at 
least (L 2 + 2L + 2)/2 different sectors. 

Now, in each ergodic sector, it is clear that the uniform state 
(equal weight superposition of all states in this sector) is a 
ground-state, that can be shown to be unique using Perron- 
Frobenius theorem. We have observed numerically on finite 
clusters that all sectors verifying 0 < |W,| + \Wy\ < L/2 
are indeed ergodic, and thus host one single ground-state 
each. However, the sectors corresponding to the boundary 
(\W X \ + \ W y \ = L/ 2 ) behave differently. As shown in Fig. 7, 
starting from a non-flippable staggered dimer configuration 
(for instance, the unique state in the sector (0,L/2)), one 
can flip dimers along an arbitrary number of diagonal lines 
to generate other non-flippable configurations (i.e. zero en¬ 
ergy states and thus ground-states) in other sectors along the 
boundary. Using simple combinatorics, one can construct 
( L p 2 ) = 0^/2) \/{p\(L/2 — p )!) configurations in the sec¬ 
tor (p, L/2 — p) for instance. 

As a result, we can construct several ground-states and we 
predict a total degeneracy of (L 2 — 2L — 6)/2 + 2 L / 2+2 . 
We have numerically checked that these are indeed the only 
ground-states (for instance 17 for L = 4, 41 for L = 6 and 
so on). In particular, the degeneracy is exponential with the 
linear size (while the number of winding sectors is only poly¬ 
nomial). This may explain why perturbing the RK point can 
lead to a variety of phenomena . 38 

This counting also provides the number of non-flippable 
states , that belong to the winding sectors along the edge, equal 
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to 4(2 l / 2 — 1). As a final comment, since we cannot define 
winding sectors for a spin model, we can nevertheless use in¬ 
formation about the momentum and/or point group symmetry 
of these states: by inspection, all these ground-states have mo¬ 
menta of the form (±fc, ±fc), i.e. along the diagonals of the 
Brillouin zone. 


1 1 

h ~ 7(0) 8 (8 - 7a 4 - a 8 ) X 
((2(Vl-a 4 -l) ~a 4 (4\/l —a 4 + 13))cos6> 

+ ^3 (4^/1 - a 4 + ll) a 4 + 6 (yi - a 4 - l)) sin 7) 

(B6) 


Appendix B: Derivation of the GQDM 


The GQDM hamiltonian is derived from the spin model us¬ 
ing the techniques described in Ref. 12. It allows us to obtain 
the values of the GQDM coefficients including all orders in 
a 2 . For conciseness we only give in this appendix the key 
ingredients for the derivation of the GQDM based on hamil¬ 
tonian Eq. (13) and refer the reader to Ref. 12 for extensive 
details on the technique. 

The first step relies in the direct evaluation of matrix ele¬ 
ments of hamiltonian Eq. (13) : 


h u = ■ cos 6 + 3 sin#) 

4n(0) 

h n = o * (cosfl + 3sinfl) 
8n(0) 

h \—1 = h 1—1 = h |—□ = hy.Q = 0 


(Bl) 

(B2) 

(B3) 


with n(0) = | cos 6*| + | sin#| 

After resummation of all contributions in a 2 , the coeffi¬ 
cients of the model are: 


1 3a 2 

** = ^W8Tr^ r ^) <cose_3i “ ,J> <B4) 

Vi = — , 1 — tt((1 + 2a 4 ) cos# + 3(1 - 2a 4 ) sin Q) 
n(6) 8 (1 — ar) 

(B5) 


1 1 

V6 ~ n(6) 8 (8 - 7a 4 - a 8 ) X 
( 5a 8 + 2 ^4\/1 — a 4 — 7^ a 4 — 4 \J 1 — a 4 + 4 ) cos 0 

+ (3 (3a 8 + (2 - 8\/l — a 4 ) a 4 - 4\J 1 - a 4 + 4 )) sin 7) 

(B7) 

, _ 1 1 

1/4 ~ W) 8 (8 - 7a 4 - a 8 ) X 

^/a 2 ^2^1 — a 4 + 7^ a 4 — \f\ — a 4 + 5 )) cos 0 

+ 3a 2 ^2 \/l — a 4 + 5 ) a 4 + \/l — a 4 — 5^ sin (9^ 

(B8) 


6 n(0) 8 (8 - 7a 4 - a 8 ) 

((3a 8 + (ll-4\/l-a 4 ) a 4 + 2( v / l-a 4 -l)) cos (9 


+ 


3a 8 + 3 (4^1 - a 4 + 3 ) a 4 + 6 (v 7 ! - a 4 - l)) sin6») 


(B9) 


Using the fact a 2 = 1/2, these coefficients take the more 
compact form: 


U = 


1 


-(cos 0 — sin#) 


4 n(6) 

V A = -7-77^ ( COS # + sin #) 
4 n(#J 


*« = ISSfe) (( 2%/5 “ 21 )+ 9 0 + 2v D sin ") 

v e = ^ ( (3 — 16X3^ cos 0 + 9 ^25 — I 6 X 3 ) sin 0 ) 

(Bl 

= 396Sm (“ X - 27 ) C “" - 9 7 - 5 ) 



(BIO) 


(Bll) 

V3^ sin#^ 


(B12) 

16 X 3 ) 

sin 0^j 


(B13) 


(B14) 


s«((“ + 8x^ cos 0 + 9 ^8x^3 - 7^ sin O^j 


792n(0) 


(B15) 
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